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■ ^ 1 INTRODUCTION 



ABSTRACT 

We describe how the large scale geonietry of the circumstellar envelopes of asymptotic 
giant branch stars can be used to probe the presence of unseen stellar companions. A 
nearby companion modifies the mass loss by gravitationally focusing the wind towards 
the orbital plane, and thereby determines the shape of the envelope at large distances 
from the star. Using available simulations, we develop a prescription for the observed 
shapes of envelopes in terms of the binary parameters, envelope orientation, and type 
of observation. The prescription provides a tool for the analysis of envelope images 
at optical, infrared, and millimetre wavelengths, which can be used to constrain the 
presence of companions in well observed cases. We illustrate this approach by exam- 
ining the possible role of binary companions in triggering the onset of axi-symmetry 
in planetary nebula formation. If interaction with the primary leads to axi-symmetry, 
the spherical halos widely seen around newly formed nebulae set limits on the com- 
panion mass. Only low mass objects may orbit close to the primary without observable 
shaping effects: they remain invisible until the interaction causes a sudden change in 
the mass loss geometry. 
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One of the striking aspects of the evolution of stars from 
the asymptotic giant branch (AGB) to the planetary nebula 
(PN) phase is the sudden change in morphology of the cir- 
cumstellar gas. On the AGB, the stellar mass loss is thought 
to be roughly spherically symmetric, but this rapidly evolves 
into prominent axi-symmetry on evolution into the proto-PN 
phase. The axial symmetry may consist of a smooth gradi- 
ent in mass loss from pole to equator or a distinct equatorial 
torus, and is often accompanied by high velocity bi-polar or 
multi-polar jets, aligned approximately along the symmetry 
axis. 

Most of the theoretical ideas that have been proposed 
to explain the axi-symmetry rely on the presence of a bi- 
nary companion. There is, however, no consensus on a spe- 
cific scenario, and exactly what happens is controversial. 
Proposals include tidal effects; engulfment of the compan- 
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ion with expulsion of a common envelope; and the effects of 
jets, from accretion disks around the companion or around 
the core of the primary (e.g., [Morris 1 987; Sokcr & Liviol 
1994 iReves-Ruiz fc Lope3ll999l : ISoker fc Rappapo7dl2000l : 
Nordhaus fc Blackmanlbood V 



If the axial symmetry of PNe is due to the effects of a 
binary companion, the relatively common occurrence of this 
phenomenon (e.g., S_ahai et al, 2007) strongly constrains the 
possible star systems that become PNe. However, the ob- 
servational situation with respect to binarity is not clear. 
Relativel y few AGB st ars are known to have binary com- 
panions |jorissenll2003l . Table 9.1). Unless the companion is 
bright or hot, it is difficult to detect given the high luminos- 
ity and variability of the primary and its thick circumstellar 
envelope. The central stars of PNe are somewhat more ac- 
cessible targets, and extensive searches have been made for 
companions using photometry, imaging, an d radial velocit y 
measurements (e.g.. ide Marco. 2006: Misz alski et al.l [20091 ). 
However, only a few dozen cases with companions have been 
found. These are typically short period systems which im- 
ply direct interaction in the past, but the general issue of 
PN formation remains an open question, especially in view 
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of the difficulty in finding companions witli intermediate or 
long periods. 

In this paper we discuss a different approach for prob- 
ing the presence of a binary companion that might affect 
the AGB-PN transition. It is based on the influence of the 
companion in shaping the mass loss of the AGB star before 
the transition takes place. The shaping effect is frozen into 
the expanding circumstellar envelope, and is reflected in its 
large scale geometry. 

This approach is possible because changes in the sepa- 
ration of the system are reasonably well prescribed. In cases 
where the companion comes into contact with or is engulfed 
by the AGB star during the transition, it must initially have 
been relatively close, within several AU. If the companion 
induces the transition from somewhat farther out, the or- 
bital separation is expected to undergo a slow increase in 
response to the mass loss of the system. In all cases where 
the companion significantly affects the transition, it is likely 
to be well within the circumstellar envelope during the AGB 
phase, and can be expected to affect the large scale geome- 
try of the envelope over a long time period. The geometry of 
the extended envelope can be studied with a variety of tech- 
niques on the AGB. It can also be studied in the halos of 
proto-PNe and PNe, in material that was ejected earlier by 
the AGB star and lies outside the developing axi-symmetric 
nebula. 

The plan of the paper is as follows. In the next section 
we discuss the effects of detached binary companions on the 
large scale geometry of AGB envelopes. We then develop 
a prescription for the observed shapes of the envelopes in 
terms of the binary parameters and the orientation. We il- 
lustrate how these results can be used to probe the presence 
of companions by measuring the shape in well-studied cases, 
and we discuss some implications for the formation of PNe. 



2 SHAPING EFFECTS OF AN AGB 
COMPANION 

There are no systematic observations of the extended en- 
velopes of AGB stars with companions, although there are 
a few cases known in which strong interactions already 
appear to be taking place. For example, tt^ Gru shows 
an equatorial torus an d jets, similar to the configuration 
seen in p roto-PNe (e.g. JChiu eral]|2006l : [Sacuto et al.ll2008l : 
iHugginsllioOTi ). We therefore rely on numerical simulations 
to underst and the effects of a companion on the large scale 
structure l|Theuns fc Jorissenlll993 : Mastrodomos & Morris 
I1999I [hereafter MM99]: lGawrvszczak et al.ll2002l '). The avail- 
able simulations are not very extensive, but the general na- 
ture of the effects caused by the companion is clear. In ad- 
dition to local interactions which may include tidal spin-up 
of the primary and accretion of envelope material by the 
secondary, there are two main effects on the extended cir- 
cumstellar envelope. 

First, for a wide range of conditions, a spiral pattern 
is imprinted on the envelope, mainly as a result of the re- 
flex motion of the primary. An example of this has been de- 
tected in the envelope of the evolved carbon sta r AFGL 3068 
l|Mauron fc Hugginj [200^ : iMorris et al] 120061 ). The spiral 
formed by this mechanism is not easy to detect. The radial 
wavelength of the pattern (A) is given by A = VP, where 



V is the expansion velocity of the envelope, and P is the 
period of the binary. Even for relatively nearby AGB stars 
the pattern is resolvable only for long periods. In the case 
of AFGL 3068 the period is ~ 800 yr and the correspond- 
ing separation ~ 100 AU, which is probably too large for 
the companion to play an important role in the AGB-PN 
transition. 

The second effect of a companion on the extended en- 
velope is to modify the mass loss geometry by gravitation- 
ally focusing matter towards the orbital plane. The latitude- 
density profile is eventually frozen into the wind flow at the 
terminal velocity in each direction, and this results in a flat- 
tened large scale structure. The evolutionary time scale on 
the AGB is much longer than the binary orbital periods of 
interest here, so the spiral pattern, if present, forms a fine 
structure on the extended envelope. Averaged over many 
spiral features, the large scale structure provides a probe 
of the presence of the companion, with a characteristic size 
scale much larger than the binary separation or the spiral 
pattern. It is this large scale structure that we quantify in 
the following sections. 



3 THE ENVELOPE SHAPE 
3.1 Characterising the shape 

We develop an approximate prescription for the large scale 
binary shaping of AGB envelopes based on the binary sim- 
ulations of MM99, which are the most comprehensive. De- 
tails of the numerical techniques, physical assumptions, and 
limitations of the s imulations are discussed at length by 
iMastrodemod (|l998l ). The models we use are listed in Ta- 
ble Al. They employ a single primary mass of 1.5 Mq with 
mass loss rates ~ 1-2 x 10~^ Mq yr~^ appropriate for AGB 
stars. The parameter space covers 0.25-2.0 Mq for the sec- 
ondary mass (Ms), 3.6-50.4 AU for the binary separation 
(d), and 5-17 km s~^ for the wind velocity at the location 
of the secondary (14). The wind from the star is assumed 
to be initially spherically symmetric. A close-in secondary 
may spin-up the primary or cause other effects which may 
enhance equatorial mass loss. In this sense the simulations 
provide a lower limit to the actual degree of axi-symmetry. 

The overall shapes of the envelopes affected by gravi- 
tational focusing vary from model to model, depending on 
the parameters of the binary. For high mass, close-in com- 
panions, the wind is highly concentrated towards the orbital 
plane, and for lower mass, more distant companions the flow 
is more spherically symmetric. At large distances from the 
star system where the wind is in steady state at the terminal 
velocity, the latitude variation of the density is frozen into 
the flow. 

We characterise the latitude variation of density in the 
extended envelope using an exponential profile. The density 
is then given by the equation: 

n(r,0) = ^exp(-0/eo), (1) 

where the r is the distance from the centre, na is the equato- 
rial density at reference radius ra, is the latitude (measured 
from the orbital plane in degrees), and 9o is the scale height. 
The exponential profile is adopted as a simple approxima- 
tion to the actual latitude variation of the five simulations 



AGB envelopes as probes of binary companions 3 



reported by MM99 in their Fig. 23 and Table 3. These pro- 
files differ in detail from model to model; the procedure used 
below matches them to an exponential at the orbital plane 
and at intermediate latitudes. For the most flattened enve- 
lope (model 1), the approximation is not an accurate fit of 
the detailed profile near the plane: the mean (absolute) devi- 
ation (averaged over latitude 0°-85°) is 0.19 dex. In this case 
however, the total variation in density from pole to equator 
is extremely large, more than a factor of 300 (2.6 dex), so 
the approximation is still a useful descriptor of the over- 
all shape. For the less flattened envelopes, which are the 
main focus of later discussion, the procedure provides quite 
good approximations to the detailed profiles. For example, 
for model 3, the exponential differs from the exact profile by 
a mean (absolute) deviation of 0.03 dex. 

To determine the value of 9o for each simulation in Ta- 
ble Al we fit equation (1) to the ratio of the density at the 
equator and intermediate latitudes 50° and 80°, and adopt 
the geometric mean. The dispersion of the individual esti- 
mates about the adopted values is 10% (rms). 

A shape parameter equivalent to 6o is the density con- 
trast (neq — npo)/npo, where Ueq and ripo are the densities 
at the equator and poles, at the same radial distance. It is 
convenient to write the contrast as Kn — 1 where Kn is the 
equator/pole density ratio, i.e., 

Kn-1 = (neq - npo)/npo. (2) 

This is zero for a spherical envelope, and increases with the 
degree of flattening. For the exponential model, it is seen 
from equation (1) that Kn and 9o are related by the expres- 
sion Kn = exp(9O/0o)- These simple characterisations are 
used to represent the overall shape of each envelope. 

3.2 Relation to binary parameters 

The next step in our prescription is to determine an empiri- 
cal relation between the binary parameters and the envelope 
shape, using the density contrast Kn — 1 introduced above. 
We follow in part the earlier discussion by MM99. 

In order to guide the analysis, we assume that the grav- 
itational focusing depends only on the following parameters: 
the masses of the primary and secondary (Mp and Ms), the 
separation (d), and the velocity of the wind at the orbit 
of the secondary (14). We ignore the possible influence of 
other quantities such as the dimensions of the stars, the 
form of the wind acceleration curve, etc., which we assume 
to be of secondary importance. It then follows from dimen- 
sional analysis that the density contrast Kn — 1 depends on 
a function $ of at most two dimensionlcss combinations of 
the independent variables and the gravitational constant G. 

In certain limiting cases, physical arguments provide 
further information on the form of For example, in 
the limit that Vo <^ Vs, where Vo is the orbital velocity 
{Va = G{Mp + Ms)/d), the trajectories of the wind parti- 
cles in the ballistic limit scale with only one dimensionlcss 
parameter, GMs/Vs d. Similarly, in the the limit of large 
Mp where the orbital velocity is large, the presence of the 
secondary acts effectively as a thin ring, and the contrast 
again depends on the same dimensionlcss parameter. In the 
intermediate regime of interest here, the situation is less 
amenable to simple analysis, but the above considerations 
motivate a parameterisation of the form: 




Figure 1. Envelope shape parameter K„ — 1 vs. binary param- 
eters Ms/V^d, in units of M0/(km s"'^)^ AU. The symbols show 
model data from Table Al and indicate the value of the param- 
eter 1 -I- Mp/Ms- open circles (1.75), filled circles (2.5), triangles 
(4.0), and squares (7.0). The solid line is the least squares fit to 
the data, and the dashed lines are one-sigma confidence limits of 
the fit. 
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Figure 2. Envelope shape parameter Kn — 1 vs. separation d, 
projected from Fig. 1 for \4 = 10 km s"'^. The upper lines and 
circles are for Mg = 1 Mq; the lower lines and crosses are for 
Ms = 0.25 Mq. 

Kn-l = HGMs/Vf d, (Mp + M,)/M,). (3) 

Ftom the numerical results of the simulations, it is found 
that K„ — l shows a strong dependence on the first parameter 
in equation (3). The data are shown in Fig. 1. The variation 
is well approximated by a power law, and a least squares fit 
gives: 
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Figure 3. Envelope shape parameter Kn — 1 vs. wind velocity 
Vs, projected from Fig. 1 for separation d = 12 AU. The upper 
lines and circles are for Ms = 1 Mq; the lower lines and crosses 
are for Ma = 0.25 Mq. 



\og{Kr^ - 1) = 5.19 (±0.66) + 1.40 (±0.21) log(A/s/V;^d).(4) 

This equation is compared to the data in Fig. 1. The solid 
line shows the fit, and the dashed lines show the one sigma 
(standard error) confidence limits. One data point (model 
12, Table Al) appears as an outlier and might be affected 
by resonance or the residual effects of a spiral shell. 

The numerical results of the simulations show no adi- 
tional variation of Kn — 1 with the second parameter of equa- 
tion (3). This is illustrated in Fig. 1 using different symbols 
to denote different values of 1 ± Alp /Ms, which range from 
1.75 to 7.0 (for = 0.25-2.0 M© and Mp = 1.5 M©). 
It can be seen that all points follow essentially the same 
curve. We can set an approximate upper limit on the influ- 
ence of this parameter by including it in the argument to the 
power law fit as GMs/K^ d x (1 ± Mp/M^)". We find that 
a = -0.29 ± 0.48 (or -0.10 ± 0.33 ff we omit the outlier) 
indicating little or no dependence on the second parameter, 
and therefore little or no dependence on the primary mass. 
Thus, even though the simulations were made for a single 
value of Mp (1.5 Mq), we can expect roughly similar results 
for other values of Mp in the relatively small mass range 
(~ 0.6 to a few Mq) relevant to stars in the later stages of 
the AGB. 

We conclude that equation (4) gives the main paramet- 
ric dependence of the envelope shape on the binary param- 
eters in the regime of interest, and we adopt this relation 
for the remainder of the paper. The dispersion of the model 
points, characterised by the median deviation, is 0.19 dex. 
This translates into a dispersion of the observed column den- 
sity contrast (Cn — 1 see below) of ~ 25%. Factors which 
may contribute to the dispersion are discussed in appendix 
A. To show the dependence of Kn — 1 on individual param- 
eters. Fig. 2 plots equation (4) as a function of d, for fixed 
values of Ms and Vs, and the model points scaled to these 



values using the equation. Figure 3 shows a similar plot as 
a function of Vs, for fixed Ms and d. Although the relations 
are not tightly constrained, it can be seen that they are suffi- 
ciently well determined to allow a preliminary investigation 
of the shaping. 

The above results are in accord with the discussion by 
MM99. They found that the quantity Ricc/d is useful for 
distinguishing different morphological types of envelopes, 
where Race is an effective accretion radius {— 2GMs/Vs) 
using the velocity of the wind at the secondary. Although 
Racc/d has a strict physical interpretation in terms of accre- 
tion only in the limit Vo <^Va, the parameterisation of Ms, 
V^, and d is the same as in equation (4) above. MM99 also 
noted that an analytic formula for the rat io of the ac c retion 
rate to the mass loss rate, called afocus bv lHan et al.l (Il995l ) 
and used by them to characterise the degree of focusing, 
is, in fact, only a fair indicator of the morphological type. 
Ofocus may be written as (GMs/d)^ /VsV^'^ , where Vr is the 
relative velocity of the wind at the location of the secondary 
(Vr^ = Vo + Vs). In the limit 14 ^ Vs, Ofocus reduces to a 
parameterisation equivalent to equation (4), but in the gen- 
eral case it differs. We find that a power law fit of Qfocus to 
K„ — 1 is a poorer predictor of the density contrast than 
equation (4), with a dispersion ~ 0.12 dex larger. Equation 
(4) is therefore the preferred relation. More extensive simu- 
lations over a larger range of parameter space are needed to 
investigate this further. 

3.3 Observed shapes 

The final step in making the complete connection between 
the binary properties and observations, is to relate the in- 
trinsic shapes of the envelopes (given above by K„ ~ 1) to 
the observed shapes. This involves a dependence on the ori- 
entation of the envelope, which we take to be the inclination 
angle (i) of the symmetry axis to the line of sight, and on 
the kind of observation. For example, it may depend on the 
distribution of column density, emission measure, or some 
other quantity, projected onto the plane of the sky. 

3.3.1 Column density 

We first consider the case where the observed intensity de- 
pends on column density. This is relevant to imaging the 
envelope in optically thin, dust scattered ambient Galactic 
light, and other techniques such as infrared imaging or mil- 
limetre molecular line observations which can be used to 
determine the column density. We characterise the observed 
shape of the envelope by the ratio between the column den- 
sity on the projected major and minor axes at the same 
distance from the centre, i.e., Cn = A'^niaj/A'^niin. It is then 
straightforward to show that for an edge-on (i — 90°) enve- 
lope: 

Ci,{90)=lnKn Kn/{Kn-l). (5) 

This relation is shown as the solid line in Fig. 4. Highly non- 
spherical envelopes with large values of K„ appear fiattened, 
with Cn(90) ~ InKn; and nearly spherical envelopes with 
values of K„ close to 1, appear nearly circular with Cn(90) ~ 
1. 

For other inclination angles, C'-N{i) has been evaluated 
numerically. The results for several values of K„ are shown 
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Figure 4. Relation between sliape parameter Kn — 1 and intensity 
contrast on the major and minor axes for envelopes seen edge-on 
(j = 90°). Solid line: Cn (intensity oc column density). Long 
dashed line: Cb (intensity oc emission measure). Short dashed 
line: Cg (scattered central illumination). 
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Figure 5. Intensity ratio of envelopes as a function of inclination 
angle. Solid lines: Cn for values of Kn = 1.5, 2.0, 4.9, and 16.8 

(bottom to top). Long dashed lines: Ce. Short dashed lines: Oq. 
For Ce and Cs, Kn = 1.5 and 4.9. 



by the solid lines in Fig. 5. In each case C^ii), given by 
equation (5) for i = 90°, decreases with decreasing inclina- 
tion and reaches 1 (i.e., the envelope appears circular) when 
the system is seen face-on. 

Wc now have a prescription for the observed shape of an 
envelope formed by gravitational focusing in a binary sys- 
tem. Given the binary parameters Ms, Vb, and d, equation 
(4) determines the intrinsic structure of the envelope, equa- 
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Figure 6. Intensity ratio Cn(90) (intensity oc column density) 
of envelopes seen edge-on, as a function of the secondary mass 
and separation for Vs = 10 km s^^. The light contours represent 
regions beyond the parameter range of the simulations. 



tion (5) gives the shape seen edge-on, and Fig. 5 shows how 
this appears at other inclinations. 

For later reference. Fig. 6 shows the variation of Cn(90) 
(the edge-on case) in the secondary mass-separation plane 
for a representative value of 14 (10 km s~^). Lighter lines 
are used to extend the contours beyond the actual range of 
the models (Table Al), assuming that equation (4) remains 
valid. The velocity of 10 km s~^ is the typical terminal ve- 
locity of an AGB wind. For small separations, the velocity at 
the secondary will be less than this (depending on the wind 
acceleration curve) and the values of Cn will be correspond- 
ingly larger; the results of specific cases can be calculated 
using the equations given. 



3.3.2 Emission measure 

If the image intensity depends on some quantity other than 
the column density, the observed shape will in general be 
different. A dependence on emission measure arises when 
the AGB halo of a newly formed PN is photo-ionised, and is 
imaged in the continuum or an emission line whose intensity 
depends on J n^dl along the line of sight. For this case it may 
be shown that Cb, the intensity ratio on the projected major 
and minor axes, is related to the shape parameter Kn by: 



Cb(90) =\nK\ 



2 J^^(l + 4(Ini^„/^)^) 
Kl - (l + 8(lnif„/7r)2)' 



(6) 



for an envelope seen edge-on (i = 90°). This relation is 
shown as the (long) dashed line in Fig. 4. It can be seen 
that imaging which depends upon emission measure pro- 
duces a higher contrast than column density, because of 
the increased sensitivity to the density. The variation of Ce 
with inclination angle has been evaluated numerically and 
is found to be similar in form to the curves for Cn- Two 
examples are shown in Fig. 5. 
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The shapes of non-spherical AGB halos around PNe 
are expected to vary over time. When the envelopes be- 
come photo-ionised, gradients in the gas pressure will pro- 
duce flows that tend to smooth out any density gradients 
l|Mellema fc FrankI [iQQSh . The smoothing time scale at ra- 
dius r is f ~ r/us where Vb is the sound speed in the gas. For 
Vs = 10 km s^^, t ranges from 3,000 to 30,000 yr for radii of 
10^'' to 10^** cm, which are the size scales typically accessible 
to observations. Thus contrast in precursor AGB envelopes 
will still be present in the AGB halos of newly formed PNe, 
but will eventually disappear in more evolved nebulae. 



3.3.3 Central illumination 

A variety of other cases can arise when the envelope is ob- 
served in dust-scattered light which originates from the cen- 
tral star or nebula. The simplest situation is optically thin, 
isotropic scattering which may be a reasonable approxima- 
tion for some halos of proto-PNe and PNe where local ab- 
sorption effects are not dominant. Other cases, in which the 
absorption of the radiation from the centre plays an im- 
portant role, need to be considered on a case-by-case basis. 
The principles, however, are the same as for those already 
discussed. 

For the optically thin scattering case, it may be shown 
that Cs, the intensity ratio on the projected major and mi- 
nor axes, is related to the shape parameter Kn for an enve- 
lope seen edge-on {i — 90°) by: 



Cs(90) =lnJf„ 



Kn{l + (lnK„/7vf) 
A'„-(l + 2(lniC„/7r)2)' 



(7) 



This relation is shown as the (short) dashed line in Fig. 4. 
It lies between the column density and emission measure 
curves. The variation of Cs with inclination angle has been 
evaluated numerically, and two examples are shown in Fig. 5. 
It can be seen that they are similar in form to the curves for 
Cn. 

For both the scattering and the emission measure cases, 
figures equivalent to Fig. 6 can be generated by replacing the 
contour values of Cn with the equivalent values of Ce and 
Cs found in Fig. 4, or directly by combining equations (6) 
and (7) with equation (4). 

In the cases where dust plays a role in the observations, 
we have implicitly assumed that the dust/gas ratio is con- 
stant throughout the envelope. In fact, there is expected to 
be a latitude variation in the dust-gas drift velocity because 
of the variation in density, and this will produce differences 
between the dust and gas envelopes. There are currently only 
limited observations bearing on the question of dust-gas cou- 
pling in AGB envelopes (e.g.. lMauron fc Huggins|[2000l ). but 
it will need to be taken into account in detailed studies. It 
should not affect that main conclusions of the present paper. 



4 DISCUSSION 

4.1 Probing the presence of companions 

The prescription given above for the large scale shaping of 
AGB envelopes by binary companions is potentially useful in 
a variety of situations because the envelopes can be observed 
with different techniques at optical, infrared, and millimetre 



wavelengths; they can also be traced in the halos of proto- 
PNe and PNe. 

In principle, the envelopes of AGB stars with known 
companions could be used to check the simulations, but in 
practice there are currently no useful observations available 
for this purpose. In the case of Mira, which first comes to 
mind, the envelope appears to be completely disturbed, from 
the outside by high velo city interaction with the interstel- 
lar medium l|Martin et al, 2007), and from the inside by the 
action of bipolar jets ( Josselin et al.l l2000l : iMeaburn et al.l 
I2009D . For other known binaries, there are no suitable obser- 
vations of the extended, undisturbed, envelopes. We there- 
fore proceed on the basis that the simulations are correct, 
bearing in mind that the prescription can be updated when 
additional information becomes available. 

One type of application is the analysis of survey obser- 
vations of the shapes of AGB envelopes in which the orien- 
tation of any binary contribution is unknown. In this case, 
population syntheses based on random orientations and as- 
sumed distributions of companion masses and separations 
will be needed to evaluate the effects on the observed shapes. 
Thus the inclination can be accounted for statistically. Since 
the median inclination for random orientations is 60° , Fig. 5 
shows that the shaping signal is well preserved for part of 
the population under these conditions. 

A second type of application is the analysis of individual 
systems in which the orientation of a possible underlying 
binary system is already suggested by observations of other 
structures such as bipolar flows. In this case the inclination 
can be directly taken into account, e.g., using the curves of 
Fig. 5. Examples of this approach are given in the sections 
that follow. 

It is interesting to note that the gravitational focusing 
expressed by equation (4) will respond to changes in the 
parameters during evolution on the AGB. For example, an 
increase in the wind velocity and an increase in separation 
due to mass loss leads to a secular decrease in the degree of 
flattening of the envelope. This effect offers diagnostic pos- 
sibilities, but is currently beyond observational capabilities. 
The sudden transition from quasi-spherical to axi-symmetry 
during PN formation (see Section 4.3), is attributable to 
other physical effects. 



4.2 IRC+10216 

As an illustration of our approach we consider the enve- 
lope of IRC4-10216, the nearest AGB star with a high 
mass loss rate. The envelope has been extensively mapped 
in th e millimetre lines o f CO and other m olecular species 
fe.g.. lHuggins et al.|[l983 : iFong et al.ll2006l ). and it has been 
imag ed at high resolution in dust scattered galac tic light 
( Ma uron fc Hug gins 1991, |2000| : iLeao et all l2006f ). These 
different observations all show that the extended envelope 
(10-200") is roughly circularly symmetric. 

In the centre of the envelope on arc-second size 
scales, HST observations reveal a bipolar structure which is 
thought to represent the early development of axi-symmetry 
which characterises t he transition to the proto -PNe phase 
(ISkinner et al. 1 Il998l : iMauron fc Huggins! |2000| ). If a com- 
panion is the cause of the bipolarity, it must have been 
present during ejection of the entire circumstellar envelope. 
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Figure 7. The envelope of IRC+10216: intensity vs. position an- 
gle for annuli of radii 10-20" (top), and 40-50" (bottom). The 
horizontal line near 165° marks a region in the data affected by 
the halo of a bright field star. The arrows indicate P.A.s orthogo- 
nal to the bipolar axis; no systematic equatorial enhancement is 
seen. 



The shape of the envelope can therefore be used to probe 
the presence of the companion. 



4-2.2 IRC+10216: limits on a companion 

The bipolar structure at the centre of IRC-f 10216 lies at 
a projected posit i on an gle of ~ 15°, and from modelling 
by ISkinner et all ()l998l ) the inclination to the line of sight 
is 70° ± 10°. The expansion time scale for this featur e is 
~ 50 yr, so it is very recent l|Mauron fc Hugging |2006| ). If 
it is caused by interaction with a companion, the compan- 
ion must have been present during the ejection of the more 
extended envelope sampled in Fig. 7. 

The inclination angle of the bipolar axis means that 
the system is seen close to edge-on. There is, however, no 
evidence for enhanced emission on the projected equator 
(indicated by the arrows in the Fig. 7). Although some of 
the arc fine-structure deviates from circularity, the average 
distribution of matter seen in the figure is relatively constant 
with position angle. The intensity ratio between the equator 
and the bipolar axis is close to 1. Even after allowing for the 
inclination, Cn l-l- 

Reference to Fig. 6 shows that these results place 
stringent constraints on any hypothesised companion to 
IRC+10216. For example, if the bi-polarity is caused by 
the onset of interaction with a companion initially within 
~ 10 AU of the primary, the possibility that it can be a 
stellar object of mass > 0.1 Mq appears to be ruled out, 
otherwise the envelope would exhibit observable large scale 
asymmetry. This result is even stronger when one considers 
a lower value of Vs appropriate for small separations, rather 
than the nominal 10 km s~^ used for Fig. 6. 

Because of the low limit of Cn, Fig. 6 also shows that 
even at distances larger than 10 AU from the primary, the 
meiss of any companion is still stringently constrained. 



4.2.1 IRC+10216: envelope shape 

We quantify the large scale shape of the envelope using ob- 
servations of dust-scattered, ambient Galactic light. Away 
from the centre, the scattering optical depth is small, and 
the imaging provides an excellent probe of the projected en- 
velope shape. The observations we use were made in the 
V band with the VLT (FORSl program 64.L-0 204, PI: de 
Laverny). The dat a were previously describe d bv lLeao et al.l 
l|2006ll , (see also iHuggins fc MauronI lioO^ ). and the pro- 
cessed image is seen in their Figs. 3 and 4. The scattered 
light images of IRC-l-10216 show substructure in the form 
of m ultiple arcs which have pr eviously been discussed in de- 
tail l|Mauron fc Hugginsll2000l ) . We are interested here in the 
large scale geometry averaged over the substructure. 

After stars and galaxies were filtered from the image, 
the envelope was divided into a series of annuli, 10" in width 
and centred on the star. The annuli were further divided into 
10° bins in position angle. The pixels in each bin were then 
averaged to determine the envelope intensity above the back- 
ground. The results are shown in Fig. 7 for two annuli which 
were the least affected by the halos of bright stars in the 
field. The formal errors in the intensity are comparable to 
the symbol size. The radial dimensions of the annuli a.re 10- 
20" and 40-50"; for a distance of 120 pc iILoup et al.lll993l') 
and an expansion velocity of 14.1 km s~^ ( Huggins fc Healvl 
1 19861 ). they sample the envelope geometry on time scales of 
400-2000 yr. 



4.3 Binary interaction in the FN transition? 

A similar situation is seen in PNe and proto-PNe where 
the axi-symmetric structure of the nebulae is often sur- 
rounded by an extended halo, which is identified with earlier 
mass loss of the star on the AGB . Numerous example s are 
seen in the s urvey observations of ICorradT et al.1 l|2003h and 
ISahai et al.l l|2007h . On the largest size scales the dominant 
shaping effect of halo gas is interaction with the interstellar 
medium (e.g., Warcing et al. _2007). Closer to the star the 
halo shape is expected to refiect the mass-loss geometry on 
the AGB. It can therefore place important constraints on 
the origins of bipolar structure in the nebulae, and we dis- 
cuss this below. It is notable that the PNe with known close 
companions lack halos (Bon d fc Livid[r99a ) . which suggests 
that these form before the late stages on the AGB. 



4.3.1 Halo shapes 

There is little detailed information on the intrinsic shapes of 
halos, but they must typically be close to spherically sym- 
metric. There are few, if any, that appear fiattened, and they 
appear round in projection even when the axi-symmetric 
structure at the centre is seen close to edge-on. The equato- 
rial tori, which often form part of the structure of PNe, do 
not extend into the halos, and there is good evidence that 
the tori develop on relativ ely short time scales along with 
the axi-symmetry and jets l|Huggins^l2007^ . 
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In the case of evolved PNe, where the nebula is com- 
pletely photo-ionized by the central star, the gas pressure 
may play a role in the evolution of the shape of the halo, 
as described in Section 3.3.2. This effect is less important 
on shorter time scales, and unimportant in proto-PNe, be- 
fore the onset of ionisation. Hence the shapes of the halos 
of proto-PNe provide direct information on the geometry of 
the mass loss before the nebula forms. 

In a r ecent HST imaging survey of proto-PNe, 
ISahai et all |20o3) found that the majority (70%) of the 
23 objects observed have detectable halos with discernible 
shapes. Essentially all of these have axi-symmetric struc- 
tures in their central regions. Most of the halos (75%) are 
observed to be round. The remainder appear elongated along 
the bipolar axes of the central structure and are most prob- 
ably illuminated by (or through) the bipolar lobes; in these 
cases the overall geometry of the halo is not revealed, but 
is consistent with spherical symmetry and non-uniform il- 
lumination. Thus the sample as a whole provides evidence 
that mass loss before the onset of axi-symmetry, is mostly 
spherically symmetric. 

In order to quantify the sh apes of the halos w e have 
examined the images reported bv lSahai et al.l (|2007| ) of sev- 
eral proto-PNe, in which the bipolar structure at the cen- 
ter has a considerable inclination angle to the line of sight, 
as evidenced by their detailed shapes and absorption lanes. 
These include IRAS 13357-6442, IRAS 17440-3310, IRAS 
18420-0512, and IRAS 19475-^3119. We characterise the 
observed shapes of the halos using the intensity ratio at 
equal distances from the center on axes aligned with the 
projected equatorial and polar axes of the central bipolar 
structure (avoiding diffraction spikes, jets, etc.). These ra- 
tios, which correspond to Cs in section 3.3.3, all lie in the 
range 0.6-1.2. The measured uncertainties are < 15%, and 
there is good agreement between the results obtained from 
images with different filters. The individual inclination an- 
gles are not known, but they are large by selection (see 
above). It is reasonable to assume that they are typically 
> 60°; using the inclination calculations of section 3.3.3 we 
can set an upper limit for the edge-on value of Cs of < 1.5. 
The corresponding column density ratio Cn is ~ 1.3. There 
is no evidence for any substantially fiattened halos lying in 
the equatorial plane of the central bipolar structures. 

The typical size scale sampled in these halo observa- 
tions is ~ 10^^ cm. For a typical halo velocity of 10 kms~^ 
the corresponding expansion age is ~ 3000 yr. The ages of 
the central bipolarity are typically less by factors of ~ 5-10 
on account of the smaller radial extents and larger expan- 
sion velocities. Such rapid time scales are consistent with 
detailed studies of other proto-PN e where the vel ocities and 
inclination angles are known (e.g.. lHugginsll2007l ). 



4-3.2 Limits on companions 

The approximate spherical symmetry of proto-PN halos and 
the sudden development of axi-symmetry in their inner re- 
gions is striking, and can be used to constrain formation 
scenarios involving binary companions. 

The scenarios which have been proposed can be di- 
vided into two general types. In the first, the formation 
of axi-symmetry and jets in the AGB-PN transition is 
the result of relatively close interaction of the primary 



with a companion, e.g., tidal spin-up by the compan- 
ion, or en gulfment follow-ed b y envelope spin-up or ejec- 
tion (e.g., 'Soker fc Livio 'l99j; iReves-Ruiz fc L6pezl Il999l : 
iNordhaus fc Blackman 2006. ) . This type of interaction can 
give rise to a n event with sh ort time scales, consistent with 
observations (|Hugginsll2007l ). In the second type of scenario, 
the components interact via accretion of th e AGB wind 
by the secondary, w hich produces jets (e.g., iMorrisI [l987l : 
ISoker fc Rappaportli200 0). In this case the jets may be trig- 
gered by enhanced mass loss from the pri mary which pro - 
duces enhanced accretion by the secondary (|Hugginsll2007l ). 

In either case, the relative sphericity of the mass loss 
before the sudden PN event severely constrains the mass 
and separation of any companion. For example, at sepa- 
rations ~ 10 AU, Fig. 6 shows that for our nominal and 
conservative limit on asphericity of 1.3, the effect of focus- 
ing alone places an upper limit of ~ 0.2 M© on the mass 
of the secondary. This conclusion is not an artifact of our 
data fitting. Models 13 and 18 (Table Al) with secondary 
masses Ms = 0.25 Mq and separations 6.3 and 12.3 AU, 
both have envelopes with large (> 3) equator-to-pole den- 
sity ratios and could readily be distinguished as non spher- 
ically symmetric. The upper limit on the mass is even more 
stringent if one considers a more realistic (lower) value for 
Vs than the nominal 10 km s~^ used for Fig. 6. In addition, 
for very close encounters, other effects are expected to en- 
hance or dominate mass loss in the equatorial plane. These 
include the infiuence of tidal effects on t he mass loss of the 
AGB star (iFrankowski fc Tvlcnda 2001), and the possible 
formation of circumbinary di sks (Van WinckeL20071 : see also 
IFrankowski fc JorissenllioOTh . 

The implications of the absence of flattened halos in 
proto-PNe and PNe seem inescapable. If relatively nearby 
binary companions are the cause of the sudden development 
of axi-symmetry and jets during PN formation, they must 
typically involve low mass stellar or sub-stellar companions. 
Only low mass objects may be present at these distances 
from the primary without the tell-tale shaping of the ex- 
tended envelope; their effect remains invisible until some 
interaction with the primary causes a sudden and major 
change in the mass loss geometry. Even at distances as large 
as 20-30 AU the absence of observable shaping of the halo 
constrains the companion mass < 0.3-0.5 M©, according to 
our nominal criterion. In fact it is an interesting and open 
question whether a companion can avoid significantly affect- 
ing the large scale shape of the AGB envelope by gravita- 
tional focusing, and at the same time accrete enough ma- 
terial to generate jets that initiate a transition to strong 
axi-symmetry. This point is currently under further investi- 
gation. 

These findings depend, of course, on the applicabil- 
ity of the underlying simulations. Those used (MM99) are 
the most detailed and comprehensive currently available. 
It is possible that effects not included in the simulations 
such as time variable mass loss (evident in the exam- 
ples discus sed above) or th e dynami cal role of magnetic 
fields fe.g.iyiemmings et al.i r2002: Hu ggins fc Manlevll2005l : 
ISabin et al. 2007 ) may affect the shaping at a quantitative 
level. Thus updated simulations exploring such effects, and 
detailed comparisons with observations of AGB stars with 
companions will be important developments for the future. 
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5 CONCLUSIONS 

This paper shows how the large scale geometry of AGB en- 
velopes can be used to probe the presence of unseen bi- 
nary companions. The presence of a companion influences 
the AGB mass loss by focusing it towards the orbital plane. 
The effect is frozen into the wind at the terminal velocity, 
and so determines the shape of the envelope at large dis- 
tances from the star. 

We have developed a prescription for the magnitude 
of the shaping effect, based on the binary simulations of 
MM99. The prescription provides a tool for the analysis of 
envelope images at optical, infrared, and millimetre wave- 
lengths. In its present form it has a number of limitations. 
Most important is the limited parameter space covered by 
the simulations. It also uses approximations to characterise 
the intrinsic shapes (Section 3.1) and in relating these to 
the binary parameters (Section 3.2). These can certainly be 
refined in the future with the help of more extensive explo- 
ration of parameter space. 

Using the prescription, we have investigated the com- 
monly observed transition from spherically symmetric mass 
loss on the AGB to axi-symmetry with jets in proto-PNe 
and PNe. If this transition is caused by interaction of the 
primary with a companion, the absence of shaping effects in 
the extended envelope (or halo) places strong constraints on 
the mass and separation of the companion. Only low mass 
objects can orbit close to the primary without observable 
effects on the envelope shape: they remain invisible until in- 
teraction with the primary causes a sudden change in the 
geometry and dynamics of ejection. 
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APPENDIX A: DETAILS OF THE MODELS 
Al Model parameters 

The simulations that we use to relate envelope shape to bi- 
nary parameters are taken from MM99. Details of the mod- 
els are given in Table Al. We omit their model 2 which 
includes the effects of spin-up, and their models 8 and 9 
which have wind velocities at the secondary that are much 
larger than expected in AGB envelopes. Table Al includes 
the binary parameters, the density ratio between the orbital 
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Table Al. Details of the binary and envelope parameters. 



Model Ms d V. V, ^ ^ eo 









(km 












01 


1.00 


3.7 


5 


10 


25.00 


200.00 


15.31 


356.63 


03 


1.00 


9.5 


11 


16 


3.20 


6.60 


42.69 


8.23 


04 


1.00 


12.6 


17 


25 


2.50 


2.90 


64.03 


4.08 


05 


1.00 


25.3 


15 


15 


1.30 


1.48 


197.20 


1.58 


06 


2.00 


37.9 


15 


17 


1.60 


1.94 


113.32 


2.21 


07 


1.00 


50.5 


15 


15 


1.50 


1.70 


136.35 


1.93 


10 


1.00 


6.3 


9 


16 


4.60 


16.00 


30.75 


18.67 


11 


1.00 


12.6 


12 


16 


1.60 


2.00 


110.81 


2.25 


12 


1.00 


12.3 


7 


10 


32.00 


83.00 


16.16 


262.13 


13 


0.25 


12.3 


7 


10 


2.60 


3.20 


59.99 


4.48 


14 


1.00 


24.6 


9 


10 


1.90 


2.70 


79.21 


3.11 


15 


0.50 


3.7 


5 


10 


7.50 


80.00 


21.28 


68.61 


16 


0.50 


9.5 


11 


15 


1.70 


2.60 


88.82 


2.75 


17 


0.50 


6.3 


9 


15 


2.90 


4.00 


52.06 


5.63 


18 


0.25 


6.3 


9 


15 


2.30 


2.60 


70.89 


3.56 



plane and latitudes 50° and 80°, and our estimates of the 
shape parameters 6o and A"„ discussed in section 3.1. 



A2 Relation of envelope shape to binary 
parameters 

The fit of equation (4) to the model data is shown in Fig. 1, 
together with the one sigma confidence limits. A higher or- 
der curve does not significantly improve the fit. There is a 
clear outlier (model 12, Table Al), so we characterise the 
scatter by the median absolute deviation, which is 0.19 dex. 
This translates into a scatter in observed quantities such 
as the column density contrast Cn — 1 (viewed edge-on) of 
~ 25%. There are several factors that may contribute to 
the scatter of the data points in Fig. 1. (1) Limitations in 
the model calculations, e.g., the influence of flnite simulation 
volumes, and different degrees to which the models reach ex- 
act steady state in the far field. (2) Variation of subordinate 
parameters, including the terminal velocity and details of 
the wind acceleration curve. (3) Additional physical effects 
such as residual spiral structure and resonance between the 
period and the wind crossing time of the orbit. (4) Devia- 
tions from the power law of equation (4), which may vary 
over the parameter space. 

The simulation grid is too sparse to explore the param- 
eter space in detail. However, it is possible to check whether 
the data are consistent with the combination of parameters 
used in equation (4), by replacing the independent variable 
with the more general form of log Ms -I- 7 log Vs + 5 log d. If we 
fit the whole data set we obtain the solution 7 = — 2.41±0.90 
and 5 = —0.50 ± 0.30, and if we omit the outlier we obtain 
7 = -2.11 ± 0.68 and 5 = -0.81 ± 0.29. The latter so- 
lution gives significantly smaller residuals and is preferred, 
but both solutions are consistent with the combination of 
parameters used in equation (3) and support the analysis of 
Section 3.2. A more extensive grid of models would be useful 
to explore the parameter space in detail. 



A3 Comparisons 

To our knowledge there are no independent simulations that 
can be directly compare d with the detaile d results of MM99. 
The models by Gawrvsz czak et al] l|2002l ) show roughly sim- 
ilar levels of gravitational focusing, but they involve some 
important physical differences concerning the geometry of 
the mass loss and the role of cooling. Unlike the mod- 
els of MM99, most of those of iGawrvszczak et al.1 (I2OO2I ) 
include initially anisotropic mass loss. The one exception 
is their model G, which has initially isotropic mass loss, 
with Mp = 1.0 M0, Ms = 0.6 Mq, Vs = 9 km s"\ and 
d — 9.6 AU. From the latitude-density profile of this model, 
wo find that Cn(90) = 1.4. This, however, is expecte d to be 
an un derestimate of the density contrast. IGawrvszczak et al.l 
(|2002l ) point out that the absence of radiation cooling in 
all their models leads to an underestimate of the density 
enhancement in the equatorial plane because the inflow is 
partly suppressed by high pressure gas that forms around 
the secondary. Consistent with this expectation, if we use the 
binary parameters of model G in equations (4) and (5), we 
predict a higher contrast (Cn = 2.4) as expected. More de- 
tailed and extensive cross-checks on the model results would 
clearly be useful. 



A4 Computational domain 

The radial extent of the computational domain of the MM99 
simulations ranges from - 1 X 10^^ to 3 X lO" cm. The 
radii of the observations of IRC-f 10216 shown in Fig. 7 are 
near the upper end of this range, and those of the proto- 
PNe halos are up to an order of magnitude larger. For en- 
velopes in steady state at the terminal velocity in each di- 
rection, the density profile remains constant with increasing 
radius, as assumed. The gas at large distances from the star 
is c ool and the turbulent v elocity (measured in few cases, 
e.g.. lHuggins fc HealvllT98^ is low, so that these should not 
have significant dispersal effects. To the extent that the sim- 
ulations have not reached exact steady state at the termi- 
nal velocity within the computational domain, there may be 
residual changes of shape farther out, but these are expected 
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to be small. Simulations on extended domains are needed to 
quantify this. 



